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We have developed a fully consistent framework for calculations in the Quasiparticle Random 
Phase Approximation (QRPA) with AW interactions from the Similarity Renormalization Group 
(SRG) and other unitary transformations of realistic interactions. The consistency of our calcula- 
tions, which use the same Hamiltonian to determine the Hartree-Fock-Bogoliubov (HFB) ground 
states and the residual interaction for QRPA, guarantees an excellent decoupling of spurious 
strength, without the need for empirical corrections. While work is under way to include SRG- 
evolved 3A r interactions, we presently account for some 3JV effects by means of a linearly density- 
dependent interaction, whose strength is adjusted to reproduce the charge radii of closed-shell nuclei 
across the whole nuclear chart. As a first application, we perform a survey of the monopole, dipole, 
and quadrupole response of the calcium isotopic chain and of the underlying single-particle spectra, 
focusing on how their properties depend on the SRG parameter A. Unrealistic spin-orbit splittings 
suggest that spin-orbit terms from the 3JV interaction are called for. Nevertheless, our general find- 
ings are comparable to results from phenomenological QRPA calculations using Skyrme or Gogny 
energy density functionals. Potentially interesting phenomena related to low-lying strength warrant 
more systematic investigations in the future. 



PACS numbers: 21.60. Jz,21.60.-n,21. 30. Fc,13.75.Cs 



I. INTRODUCTION 

In recent years, rare isotope beams have become a 
major focus of the experimental nuclear physics com- 
munity. Using such beams, more and more exotic nu- 
clei become accessible experimentally that exhibit novel 
structural features and excitations modes like neutron- 
skin vibrations and allow studies of sensitive details of 
the nuclear interactions and the theoretical models that 
are used to describe them. Over the past few decades, nu- 
clear structure calculations for medium- and heavy-mass 
nuclei were almost exclusively carried out in the frame- 
work of Density Functional Theory (DFT), using phe- 
nomenological energy density functionals (EDFs) of the 
Skyrme or Gogny type [1], or in Relativistic Mean-Field 
Theory [2] . While phenomenological EDFs uniformly de- 
scribe the bulk properties of nuclei near the valley of sta- 
bility very well, there is a significant model dependence 
and deterioration of quality for results in exotic nuclei 
and spectroscopic observables in general. While DFT 
formally resembles the Hartree-Fock (HF) and Hartree- 
Fock-Bogoliubov (HFB) methods, the EDF parametriza- 
tion contain correlation effects beyond the mean field; on 
the one hand, this allows for a better description of exper- 
imental data in a comparatively simple framework, but 
on the other hand, there is no clear connection to the un- 
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derlying AW (and 3N, 4N, . . . ) interactions, and, there- 
fore, no way to improve the EDFs in a systematic fashion, 
e.g., by many-body perturbation theory (MBPT). 

One way to overcome these problems is to use realis- 
tic AW interactions like Argonne VI 8 [3] or the poten- 
tials from next-to-next-to- next-to-leading order (N3L0) 
of chiral effective field theory (EFT) [4, 5], which accu- 
rately describe AW scattering data. The latter are par- 
ticularly appealing because chiral EFT provides a con- 
sistent set of accompanying 3A~ interactions, although 
the 3N interaction has thus far only been derived to or- 
der N2L0 [4]. Since realistic interactions induce strong 
short-range correlations in the AW system, one needs to 
tame their short-range behavior, preferably by means of 
a unitary transformation which automatically preserves 
the AW observables. Examples of such unitary trans- 
formation techniques are the Unitary Correlation Op- 
erator Method (UCOM) [6], and the Similarity Renor- 
malization Group (SRG) approach [7] , which will be the 
method explored in the following. The SRG evolution 
drives the two-body interaction to band-diagonality in 
momentum space, thereby decoupling low and high mo- 
menta. This decoupling results in soft interactions with 
greatly improved convergence properties in (quasi-)exact 
many-body methods. In addition, SRG-evolved interac- 
tions yield bound nuclei already at the mean-field level, 
and are suitable for low-order MBPT treatments (see Ref. 
[7] and references therein). 

Soft AW interactions have been discussed in nuclear 
theory for many decades, but were originally discarded 
due to their inability to produce the proper saturation 
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behavior in nuclear matter [8]. From the modern point 
of view, this merely shows that for each JVJV interaction, 
consistent 37V forces are required to properly describe 
nuclear systems. According to general EFT principles, 
up to A-body interactions must be considered in the A- 
nucleon system, for which chiral EFT guarantees a natu- 
ral hierarchy with 7V7V > 37V > 47V etc.: the leading 3TV 
force appears at N2LO, the leading 4TV force at N3LO, 
and so on. The SRG evolution beautifully illustrates the 
inseparability of the nuclear interactions, because many- 
nucleon forces are naturally induced during the SRG flow 
[7]. In the case of SRG-evolved TV TV interactions, the in- 
clusion of repulsive 37V interactions is essential to prevent 
overbinding in heavier nuclei. Formally, the 37V interac- 
tion must then be evolved consistently along with the 
7V7V interaction, which has recently been accomplished 
by Jurgenson et al. [9]. 

With realistic interactions and similarity transforma- 
tion techniques, significant progress has been made to- 
ward a comprehensive description of nuclear structure 
all across the nuclear chart. Ab initio calculations for 
light nuclei (see [6, 7] and references therein) are comple- 
mented by mean-field based approaches in heavy nuclei. 
In recent years, we have developed a framework for using 
effective interactions, given in terms of their harmonic- 
oscillator matrix elements, in a wide-range of mean-field 
based approaches, from HF and HFB to the Random 
Phase Approximation (RPA) and its extensions [6, 10- 
12]. 

The purpose of this article is twofold. First, we ex- 
tend the description of excitations to open-shell nuclei 
by means of a fully consistent Quasiparticle RPA [13, 14] : 
The QRPA is built on the ground states obtained from 
the HFB method, and the same intrinsic Hamiltonian, 
including the exact Coulomb interaction, is used in both 
the HFB and QRPA calculations. Second, we perform 
a survey of the response of isotopic chains using SRG- 
evolved TV TV interactions. While work is under way to 
include SRG-evolved chiral 3TV interactions in our over- 
all framework, a 3TV contact force (or equivalent density- 
dependent interaction) is implemented as an intermedi- 
ate step, allowing us to carry out preparatory studies and 
identify issues in anticipation of the full 3TV interaction. 
Our nuclear structure results based on SRG-evolved TVTV 
(and eventually 3TV) interactions will provide important 
guidance for ab initio DFT efforts in the framework of the 
Universal Nuclear Energy Density Functional (UNEDF) 
project [15]. 

This article is organized as follows. In Sec. II we re- 
view the QRPA formalism and provide details and tests 
of our implementation. In Sec. Ill, we fix the free pa- 
rameter of the density-dependent interaction and discuss 
some open issues pertaining to 37V forces and proceed to 
summarize HFB results for the calcium isotopic chain in 
Sec. IV. This sets the stage for the QRPA results, which 
are presented in Sec. V. Explicit expressions for the 
QRPA matrix elements are collected in the appendices. 



II. QRPA FORMALISM AND 
IMPLEMENTATION 

A. Quasi-Particle Random Phase Approximation 

Our starting point is the intrinsic Hamiltonian 

H = T int + V, (1) 
where the intrinsic kinetic energy is defined as [13, 16] 

r '»=r-T m =(l-£)££L-^£p,p M 2) 

We formulate the QRPA in the canonical basis of the 
Hartree-Fock Bogoliubov (HFB) ground state [13, 17, 18]. 
Normal-ordering the Hamiltonian w.r.t. the HFB vac- 
uum, we obtain 

H = E + J2Hlkala k ,+V ies (3) 

kk' 

where Eq is the energy expectation value in the HFB 
vacuum, and {0^,0^} are quasiparticle operators in the 
canonical basis. The residual interaction is given by 

V Ies = J ^2 Vkk'W ■ c\c\,C V Ci \ , (4) 
kk'W 

where cj. are the creation operators of the canonical basis 
in particle representation, and Vkk'w denotes an antisym- 
metrized but not normalized two-body matrix element. 
Details on how to obtain the quasiparticle representation 
of Vi-es can be found in Refs. [13, 14]. 

Assuming spherical symmetry, the canonical basis 
states come in pairs { m M )} which are related 

by time reversal: 

\Wn) = {-l) l+3 - m Vi-rn), (5) 

where fi = (nljr) is a collective index for the radial, an- 
gular momentum, and isospin quantum numbers. In the 
canonical basis, the Bogoliubov transformation between 
the particle and quasiparticle representation reduces to 
the BCS-likc form [19, 20] 

a lim = U t^ C lim + v (J, c fJ.m j (6a) 

where we have expressed the annihilation operators as 
spherical tensors [21], 

= (-iy +m a^ m = -{-lYctpm, (7) 

and absorbed a factor (—1)' into the coefficients to 
simplify the formulas. 

The QRPA phonon creation operator in the canonical 
basis has the general form [22] 

0\ — \ " X k , ,ryt ryt —Y k , ,a , ,n 

/ j y ^fj,m,ii'm' ul fj,m ul fj,'m' f+m^'m 1 U C m 'w-pm j 
(fj,m)<(fj,'m') 

(8) 
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where the sum over quasiparticle states must be re- 
stricted to avoid double counting. Since we assume spher- 
ical symmetry, it is convenient to switch to an angular- 
momentum coupled representation [14]: 

®kJM = ^im'^-nn'JM — Ypjl'-AppijM i (9) 

where the coupled quasiparticle-pair creation operator is 
defined as 



A 1 

■^■HH'JM — 



1 



W m,m 



ti m f m '\ JM ) a tm a l> m > 



(10) 

and A^ijm is its spherical adjoint [cf. Eq. (7)]. 

Using the Equations-of-Motion method [14, 22], one 
can define the QRPA matrices A and B via the commu- 
tators (fx < fj,', v < u') 

= <*| [A^jm, [H,Al, JU }} |*) , (11a) 



B 



JM — 



*l [A^jm,[H,A^,j M }} I*) 



(lib) 



where we resort to the usual quasi-boson approximation 
by assuming that the many-body state \^) is the HFB 
vacuum. For spherically symmetric systems, the QRPA 
matrices and the amplitudes X and Y are independent of 
the angular momentum projection, and one obtains the 
following reduced set of QRPA equations: 



A J B j 

-B j * -A u 



X kj 



(12) 



where hu>k is the excitation energy of the fcth QRPA state 
w.r.t. the ground state. Explicit expressions for the ma- 
trices A and B can be found in Appendix A. 



B. Transition Operators 

For electric multipole transitions, the reduced transi- 
tion probabilities arc defined as 



1 



B(EJ, Ji -+ J f ) eee ^— - \(fj f \\ Qj \\iJi ) 



(13) 



In the QRPA, we consider transitions from the + ground 
state of an even-even nucleus to an excited state de- 
scribed by the QRPA phonon operator (9), and the re- 
duced matrix element can be evaluated to (see, e.g., Ref. 
N) 



<fcj||Qj||0) 

= J2 / 1 + s ( U *»V + (-^'w) 



x (x%r(fi\\Qj\\»')H-iyY$*(vL\\Qj\\vL')*) . 

(14) 



In the limit of small momentum transfer, the multipole 
transition operator is defined as the sum of the isoscalar 
and isovector operators 



>-Y,r(Yj M (h) 



and 



Q I j V M = hT,^r{Y JM {f i ) 



(15) 



(16) 



Exceptions are the monopole operator, which would be a 
constant unable to cause transitions, and therefore needs 
to be defined as 

A 

Ooo = 1(1 + rf)r^oo(n) , (17) 



and the isoscalar and isovector dipolc operators, which 
are corrected for center-of-mass effects [13]: 



Qim — 



A 

i=l 



Yi M (fi 



and 



n iv 



N Z Z N 

e— r p Y 1M (fp) - e— r n Y 1M (r n ) ■ 



(18) 



(19) 



p=i 



n=l 



where i? ms is the intrinsic mean-square radius operator 
[23]. 

Since the NN interactions used in this work are ob- 
tained by means of an SRG evolution, we need to ad- 
dress the issue of evolving observables in a consistent 
fashion. However, it has been demonstrated in related 
approaches like the UCOM [10] and the Lee-Suzuki trans- 
formation in the No-Core Shell Model [24] that the ab- 
solute values of the transition operator matrix elements 
entering (13) hardly change. The reason for this is the 
long-range, low-momentum character of the r J operator, 
while the SRG, UCOM, and Lee-Suzuki transformation 
modify the short-range, high-momentum matrix elements 
[25]. Since there are more significant uncertainties due to 
the A-dependence of the interaction and the simple na- 
ture of the phenomenological 37V interaction, we content 
ourselves with using the uncvolved transition operators 
in the following. 



C. Calculation Details 

Our QRPA implementation is suitable for use with 
NN Hamiltonians with or without a density-dependent 
NN (or contact 3N) term. In the present study we 
use an SRC-evolved Argonne V18 interaction supple- 
mented with a phenomenological density-dependent two- 
body term to take into account missing genuine and in- 
duced thrcc-nuclcon interactions. The matrix elements 
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FIG. 1. (Color online) Number operator response for nonspu- 
rious monopole states in 56 Ca (see text). (Vsrg+DDI with 
A = 2.02 fm~\C 3 jv = 3.87 GeVfm 6 .) 



of this potential, which wc will denote as Vsrg j an d the 
two-body part of the intrinsic kinetic energy (2) are eval- 
uated in a relative spherical HO basis and transformed 
to the jj-coupled single-particle basis by using Talmi- 
Moshinsky brackets, as described in Rcf. [26]. The ma- 
trix elements of the density-dependent interaction, which 
is described in Sec. Ill, can be evaluated directly in the 
jj-coupled basis (see, e.g., Ref. [27] and also Appendix 
B). 

The input ground states for our QRPA calculations 
are obtained using the spherical HFB implementation de- 
scribed in Rcf. [27]: single-particle states are expanded 
in a spherical HO basis of 15 major oscillator shells, en- 
suring converged ground-state energies. In calcium iso- 
topes, the typical ground-state energy gains from increas- 
ing the size of the single- particle basis from 11 to 13 ma- 
jor shells are about 200 keV, and from 13 to 15 major 
shells 100 keV or less, corresponding to less than 0.1% 
of the total energy. In the tin region, the absolute gains 
arc roughly twice as large, but the relative accuracy is 
similar. The ground-state energies are minimized w.r.t. 
the oscillator length auo by considering a mesh of values 
ranging from 1.5 to 2.40 fm with a spacing of 0.05 fm. 
For 15 oscillator shells, the ground-state energies of the 
adjacent ano mesh points differ from those at the min- 
ima by 5 keV or less in the calcium chain, and 40 keV or 
less in the tin chain. 

After constructing the canonical basis of the HFB 
ground state, we determine all possible two-quasiparticlc 
(2qp) configurations for a given J 7r . We stress that the 
2qp basis is not truncated in any way, which leads to an 
excellent decoupling of spurious states in our calculations 
(see Sec. II D). To estimate uncertainties due to the dis- 
crete oho mesh, we perform QRPA calculations for the 
HFB solutions at the neighboring points as well. The 
variation of both individual excited state and centroid 
energies is about 50 keV or less in the results presented 
in the following. 
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FIG. 2. (Color online) Isoscalar (left) and isovector (right) 
dipole strength distributions in 56 Ca with ( ) and with- 
out ( ) c.o.m. correction. The discrete strength dis- 
tributions have been folded with a Lorentzian of width V — 
1.0 MeV. In the isoscalar channel, the spurious state has been 
removed explicitly from the response of the uncorrected op- 
erator (see text). (1/ S rg+DDI with A = 2.02 fin -1 , C 3 iv = 
3.87 GeVfm 6 .) 



D. Spurious States 



In a QRPA calculation, spurious states emerge as a 
consequence of the breaking of the symmetries of the 
nuclear Hamiltonian by the ground-state wavefunction. 
In the + channel, nonvanishing neutron and/or proton 
pairing break the U(l) symmetries associated with neu- 
tron and proton number (or, alternatively, nucleon num- 
ber and charge) conservation. The use of the same inter- 
action in the particle-hole and particle-particle channels 
in both the HFB and QRPA calculation ensures that the 
corresponding spurious state(s) are well-decoupled from 
the excitation spectrum [17, 18]. 

As a typical example, we show the number operator 
response for the nonspurious + states of 56 Ca in Fig. 1, 
which is less than 10 -4 / MeV in magnitude overall, and 
less than 10 -6 / MeV in the energy range up to 40 MeV, 
which is relevant for the giant monopole resonance. The 
number operator response vanishes to machine accuracy 
if the pairing collapses (i.e., in the HF + RPA limit). 
The spurious + state in 56 Ca is found at 0.1 keV, and 
the largest energies we found in our calculations are ap- 
proximately 20 keV, which still indicates excellent con- 
sistency. 

The decoupling is realized just as well for the spurious 
1~ state associated with the breaking of translational in- 
variance. To verify that our QRPA solutions are free of 
ccnter-of-mass contamination, wc compare the isoscalar 
1~ strength distribution obtained with (18) and the un- 
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corrected operator 

A 

Q{ S M = e^r!Y 1M (ri)- (20) 

i 

If translational invariance is properly restored for the 
QRPA solutions, the strength distributions must agree 
for the nonspurious states, so the only effect of the cor- 
rection term in (18) is the removal of the spurious state 
associated with the translation of the whole nucleus. 

In the isovector case, the corrected dipole operator (19) 
is used, which is equivalent to 

Z Z N 

D'-^e(r p -fi) = -e^r p --e^r„. (21) 

p p n 

If translational symmetry is properly restored, the c.o.m. 
operator R cannot cause transitions, and thus for intrin- 
sic excitations the matrix elements of D' must be iden- 
tical to those of 

z 

r> = 5>rv (22) 
p 

Figure 2 demonstrates that isoscalar and isovector 
strength distributions of the corrected operators are prac- 
tically identical to those of the operators (20) and (22), 
respectively. This confirms our previous findings for 
closed-shell nuclei [10] in the more general HFB+QRPA 
framework for open-shell nuclei. 

As in the monopole case, the spurious 1~ state itself 
lies at very low energies, independent of the mass. It 
is found between 1 and 5 keV in closed-shell nuclei, and 
below 20 keV for open-shell nuclei in all cases. These en- 
ergies are considerably lower than the spurious 1~ ener- 
gies of several hundred keV which are reported for other 
consistent QRPA approaches in the literature [18, 28- 
30]. The reason is the use of the intrinsic kinetic en- 
ergy (2) in our HFB+QRPA calculations. If we do not 
subtract T cm , the spurious state energies increase to the 
sizes reported by other groups, but the quality of the 
translational-symmetry restoration is not affected. 



III. THE DENSITY-DEPENDENT 
INTERACTION 

In the present study, we supplement an SRG-evolved 
Argonnc V18 NN interaction by a phenomcnological 
density-dependent two-body term to account for miss- 
ing genuine and induced 3iV interactions. As discussed 
in Ref. [31], the linearly density-dependent interaction 
(DDI) 

vM = ^(l + P,)p{^)s S (ri-r 2 ), (23) 





Affm" 1 ] 


C* 3 iv(A)[GeVfm 6 ] 




1.78 


4.41 


Ca 


2.02 


3.87 




2.40 


2.94 




1.78 


4.95 


Su 


2.02 


4.35 




2.40 


3.42 



TABLE I. Running coupling strength C3jv(A) for various val- 
ues A used in this work (see text). 

where P a is the spin-exchange operator, gives the same 
contribution to the ground-state energy as the 37V con- 
tact interaction 

v 3 = C 3N 5 3 (n - r 2 ) S 3 (r 2 - r 3 ) (24) 

in systems with time-reversal invariance. Such a contact 
force was recently used in conjunction with similarity- 
transformed interactions in Hartree-Fock and MBPT cal- 
culations [32]. 

Since the QRPA is the limit of low-amplitude motion 
of time-dependent HFB, the QRPA matrix in Eq. (12) 
is (up to the metric) the stability matrix of the ground- 
state energy functional [13], and the replacement of the 
3N contact term with the DDI (23) is meaningful in 
this context, provided one properly takes rearrangement 
terms due to the density-dependence into account (see 
Ref. [33] and Appendix B). 

In the present work, we treat the strength of the 
density-dependent term as a running coupling constant 
which depends on the SRG parameter A, since it is an ef- 
fective parametrization of initial 3N interactions, which 
should be present in the "bare" nuclear Hamiltonian, as 
well as induced 3iV interactions, which arc generated dur- 
ing the SRG flow (see Ref. [7] and references therein). 
We fix C3at(A) in Hartree-Fock calculations by fitting the 
experimental charge radii of a set of closed-shell nuclei, 
because the radii are much less sensitive to many-body 
corrections than the ground-state energy. 

In Table I, we list the C3jv(A) values for three A's used 
in the following; Fig. 3 shows the corresponding ground- 
state energies (per nucleon) and charge radii for the fit 
nuclei. We find that for these nuclei, the A-dependence 
of the charge radii can be absorbed into C3at(A) at the 
HF(B) level, while a many-body approach beyond the 
mean-field and a complete treatment of the 37V interac- 
tion is required to reduce or remove the A-dependence of 
the ground-state energies. 

The strength of our DDI is notably lower than the 
strength of the density-dependent term of past and cur- 
rent Skyrme functionals (see Ref. [36] and references 
therein). In part, such a strong repulsive DD term is re- 
quired to counteract the attractive terms in the Skyrme 
functional that are necessary to fit binding energies and 
radii at the same time in a mean-field calculation, which 
differs from our fit strategy 
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FIG. 3. (Color online) Ground-state energy per nucleon 
and charge radii of closed shell nuclei for Vsrg+DDI with 
(Affm- 1 ], C 3 n [GeVfm 6 ])= (2.40, 2.94) (•), (2.02, 3.87) (■), 
and (1.78, 4.41) (♦). Experimental values [34, 35] are indi- 
cated by black bars. 



We also note that we do not need to use fractional den- 
sity dependencies for the DDI to obtain reasonable GMR 
energies in Sec. II A. Modern Skyrmc functionals use 
fractional density dependencies because the strong linear 
density-dependent terms tend to overestimate the GMR 
energy, i.e., the incompressibility of nuclear matter. Since 
phenomenological Skyrme functionals with constant co- 
efficients correspond to nuclear contact interactions, non- 
integer density dependencies parametrize physical effects 
from finite-range NN interactions in addition to the 3N 
contact term [37]. 

Inspecting Fig. 3, we note that the quality of the fit to 
the charge radii varies with A: the charge radius of 16 is 
slightly too large, while the charge radii of heavier nuclei 
are underestimated. This is evidence that the limited 
spin-isospin dependence of the DDI (23), which only acts 
in the (S, T) = (1, 0) channel, is insufficient. Explicit spin 
and isospin degrees of freedom in the 3iV interaction — 
as in the chiral 3iV interaction, for instance — would lead 
to non-vanishing DDI matrix elements in the other (S, T) 
channels [38], and allow for a different scaling with the 
particle number. Since we will be focusing on the calcium 
isotopic chain and selected tin isotopes in the remainder 
of this work, we list two sets of DDI parameters in Table 
I that were optimized for the corresponding mass regions 
by fitting the charge radii of 40 Ca and 114 Sn, respectively. 

A serious issue emerges for the spin-orbit splittings, 
which are listed for various nuclei in Table II. Us- 
ing just a two-body SRG-evolved NN interaction with 
A = 2.40 fm , the spin-orbit splittings are about 2-3 
times as large as values extracted from experiment. HF 
calculations for Vsrg with A = 1.78 and 2.02 fnU 1 col- 
lapse and are not included. The repulsive DDI stabilizes 
the HF(B) calculations and leads to a compression of 
the single- particle spectra, but it results in a significant 



Level 




(A[fnT 


C 3 iv [GeVfm 6 ]) 




Exp. 




(2-4, -) 


(2.4, 2.94) (2.02, 3.87) (1.78, 4.41) 




16 ttOp 


9.61 


3.15 


3.12 


3.02 


6.32 


uOp 


10.04 


3.18 


3.14 


3.05 


6.18 


40 Ca iv Od 


15.19 


4.21 


4.14 


4.04 


6.00 


vrO/ 


12.07 


3.78 


4.10 


4.18 


4.95 


uOd 


15.95 


4.27 


4.19 


4.08 


6.00 


1/0/ 


14.90 


4.39 


4.54 


4.53 


4.88 


48 Ca vOf 


15.96 


3.51 


3.46 


3.36 


8.97 


100 Sn7rO S 


16.46 


2.45 


2.43 


2.39 


6.82 


vOg 


17.63 


2.03 


1.99 


1.96 


7.00 


132 Sn7rO# 


15.61 


1.82 


1.71 


1.67 


6.08 


uOh 


20.12 


2.98 


2.89 


2.48 


6.53 



TABLE II. Proton (tt) and neutron (v) spin-orbit splittings in 
MeV for different interaction parameters (A,C3jv), compared 
to experimental values [39, 40]. 



underestimation of the spin-orbit splittings. As a conse- 
quence, there is a shift in the major shell closures: 48 Ca, 
for instance, has a small non-vanishing neutron pairing 
energy because the O/7/2 level lies close to the pf ma- 
jor shell, while 120 Sn is essentially a closed-shell nucleus 
because the Ohn/2 level is not shifted to sufficiently low 
energies to produce the N = 82 major shell closure. 

The resulting proximity of levels with Aj = Al = 2 in 
the region of the Fermi surface favors strong quadrupole 
interactions, and causes the nucleus to develop a defor- 
mation. In our spherical QRPA calculations, this insta- 
bility w.r.t. to quadrupole deformations is signaled by 
the lowest 2 + energy becoming purely imaginary, and we 
have to discard results based on the unstable spherical 
ground-state configuration. The calcium isotopic chain is 
free of this pathology, hence we focus on these isotopes in 
the following, and defer calculations for heavier isotopic 
chains to a future publication. Since the present work 
sets the stage for the use of the chiral 3iV interaction (or 
a density-dependent variant), which contains additional 
spin-orbit and tensor terms (see, e.g., Refs. [38, 41, 42]), 
we expect a significant impact on the spin-orbit physics, 
but it remains to be seen whether the discussed problem 
can be resolved. 



IV. GROUND-STATE PROPERTIES 

Before we present results from our QRPA calculations, 
we briefly discuss the properties of the calcium ground 
states on which the QRPA is built. All results have been 
obtained with SRG-evolved Argonnc V18 interactions, 
supplemented by the adjusted DDI discussed in Sec. III. 
Odd nuclei have been treated in a self-consistent Equal 
Filling Approximation (EFA) [43]. 

In Fig. 4, we summarize the HFB ground-state proper- 
ties of the studied calcium isotopes. While the climina- 
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FIG. 4. (Color online) Ground-state properties of the calcium isotopes for Vsrg+DDI with (A[fm _1 ], C3jv[GeVfm 6 ])= (2.40, 
2.94) (•), (2.02, 3.87) (■), and (1.78, 4.41) (♦): (a) ground-state energies per nucleon, (b) charge radii, (c) odd-even mass 
differences, and (d) pairing energies. Experimental values [34, 35] are indicated by black bars or crosses. 



tion of short-range correlations via the SRG evolution 
yields bound nuclei for Vsrg+DDI at the HFB level, 
ground-state energies shown in Fig. 4a are underesti- 
mated by 4-6 MeV per nucleon, depending on the cho- 
sen A. As discussed, e.g., in Rcf. [6], this is due to 
long-range correlations that are not taken into account 
by HFB, but can be recovered to a large extent by low- 
order MBPT due to the perturbative character of the 
SRG-evolvcd NN interaction [7]. 

Figure 4b displays the charge radii, which are reason- 
ably close to experimental values due to the fit of the DDI 
discussed in Sec. III. Overall, our calculations fail to re- 
produce the pattern of the experimental charge radii be- 
tween Ca and 48 Ca, whose proper description requires 
the inclusion of effects beyond the mean- field (see, e.g., 
Ref. [44] and references therein). For 40 Ca and 48 Ca, 
though, they are not far off. 

In Figs. 4c and 4d, we show the experimental odd-even 
binding energy differences 

(if 

A (3) = 1 ; (E{N+l)-2E(N) + E{N-l)) (25) 

and the pairing energies (expressed in the canonical basis) 
[13] 

iW = ^ H( 2 ^ + l)A MM u^ , (26) 

respectively. A' 3 ), in particular, is sensitive to the single- 
particle structure of the HFB vacuum state near the 
Fermi surface. 

We find that the agreement with experimental values 
is reasonable in the range 16 < N < 28, although the 
staggering is less pronounced in our calculation, which is 
at least in part due to the lack of time-reversal symmetry 



breaking in our EFA treatment of odd nuclei. The sharp 
jump in the A^ 3 ) at N = 20 is a signal of the major shell 
closure, which leads to the expected collapse of HFB pair- 
ing in 40 Ca. A similar jump is expected at the N = 28 
shell closure, but here the experimental data are not re- 
produced, and we note that the pairing energy does not 
vanish in 48 Ca. A look at the canonical single- neutron 
energies shown in Fig. 5 helps to clarify this issue. In 
our calculation, the spin-orbit splitting of the 0/ levels 
in 48 Ca is about half the experimental splitting. Con- 
sequently, the O/7/2 level lies close to the main p/-shell, 
and the density of unoccupied levels is sufficiently high to 
produce pairing. We also note that the underestimation 
of the spin-orbit splittings is strongly isospin dependent, 




theo. exp. theo. exp. 



FIG. 5. (Color online) Canonical neutron single-particle 
energies in 40 Ca and 48 Ca for Vsrg+DDI with A = 
2.02 fm _1 , Csn = 3.87GeVfm 6 . Experimental values have 
been taken from Refs. [39, 40]. 
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FIG. 6. (Color online) Canonical single-neutron energies (left) and isoscalar monopole response (right) of 120 Sn for Vsrg+DDI 

with (A[fm" 1 ],C 3 iv[GeVfm 6 ])= (2.40, 3.42) ( ), (2.02, 4.35) ( ), and (1.78, 4.95) ( ). Results obtained with 

Gogny D1S are included for comparison ( ), and the arrow indicates the experimental centroid energy [45]. The discrete 

ISM strength distributions have been folded with a Lorentzian of width V — 1 MeV. 



because the splittings in he N = Z nucleus Ca are re- 
produced fairly well by comparison (also cf. Table II). 
As indicated in Sec. Ill, the issue of missing spin-orbit 
strength is not limited to the calcium isotopic chain or 
specific sets of interaction parameters A and C^n, and 
it becomes an obstacle for QRPA calculations in heavier 
nuclei. 

We conclude this section by mentioning that the 
present results for the odd-even binding energy differ- 
ences are compatible with those of "hybrid" studies where 
Vsrg or Mow-k are used as pairing interactions in con- 
junction with phenomenological energy density function- 
als like SLy4 [46] or Gogny D1S [27]. The inclusion of the 
DDI overcomes the low level density that obstructs pair- 
ing in HFB calculations when only similarity-transformed 
two-body interactions are used [27], and brings the single- 
particle energies closer to EDF results. Since the DDI 
(23) vanishes in the pairing channel by construction, the 
A' 3 ' and pairing energies shown in Figs. 4c and 4d probe 
only the pairing properties of Vsrg, an d we obtain mag- 
nitudes that are similar to the hybrid approaches, aside 
from issues with the systematics due to the missing spin- 
orbit strength. It will be interesting to see how the situ- 
ation changes once a consistent set of SRG-evolved NN 
and 3N interactions is used in the calculation, because 
the latter will directly affect the pairing channel. 



V. NUCLEAR RESPONSE 

A. General effects of the SRG Evolution 

We first discuss the general effects of the SRG evolu- 
tion on the nuclear response, using the isoscalar + re- 



sponse of Sn as an example. Sn is chosen because 
it has a richer single-particle structure around the Fermi 
surface than the calcium isotopes, which facilitates the 
discussion, and because its ground state is stable against 
quadrupole deformation in our calculations (cf. Sec. III). 

In Fig. 6, we show the canonical single-neutron en- 
ergies and the isoscalar monopole response function for 
Vsrg+DDI with the three parameter sets that have been 
optimized for the tin isotopic chain (cf. Table I). We 
note that the shell structure of the single-neutron ener- 
gies becomes more pronounced as A is lowered. There 
are two distinct effects: During the SRG evolution, sub- 
shells move closer together, while the gaps between ma- 
jor shells increase as the SRG evolution renders the in- 
teraction increasingly nonlocal and thereby reduces the 
effective mass [6, 23, 27]. For comparison, we have also 
calculated the single-neutron spectrum of 120 Sn using the 
Gogny D1S functional. Except for the underestimation of 
the 0<?, Oh, and Oi spin-orbit splittings, the single-neutron 
spectra of the three Vsrg+DDI and D1S are rather sim- 
ilar for the levels around the Fermi surface. The similar 
level density for Vsrg+DDI and for Gogny D1S, which 
corresponds to an effective mass of about 0.7 times the 
bare nucleon mass, suggests that the inclusion of the re- 
pulsive DDI compensates for the extremely low effective 
mass that is a common feature of low-momentum inter- 
actions. 

The isoscalar monopole (ISM) response functions dis- 
played in Fig. 6 reflect the evolution of the single-particle 
spectra with the A. The significant fragmentation of the 
ISM strength distribution for A = 2.40 fm _1 is reduced 
notably due to the bunching of levels within the major 
shells as A is lowered. At the same time, the gap be- 
tween major shells increases, leading to a shift of the 
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FIG. 7. (Color online) Running energy-weighted sums of 120 Sn for Vsrg+DDI with (A[ fin -1 ], C 3 n [ GeV fm 6 ])= (2.40, 3.42) 

( ), (2.02, 4.35) ( ), and (1.78, 4.95) ( ). The light gray dashed lines indicate the values of the classical sum rules 

(see text). 



strength distribution to higher energies. This shift to 
higher energies competes with the attractive isoscalar 
residual interaction, which aims to shift the response to 
lower energies. To get some insight into how the two ef- 
fects depend on A, we list the centroid energies for the 
unperturbed HF response without the residual interac- 
tion, and the QRPA response in Table III. The inte- 
gration intervals for the moments m /i have been cho- 
sen by inspecting plots of the response functions. In the 
HF case (not shown), the bulk of the ISM strength lies 
between 10 and 30 MeV, with a tail containing a few 
percent of the energy-weighted sum rule extending up to 
roughly 45 MeV. In the full QRPA calculation, the ISM 
strength is concentrated in a much smaller range from 10 
to ~ 25 MeV due to the attractive residual interaction 
(see Fig. 6), and the centroid lies 8.3 to 8.4 MeV be- 
low the centroid of the unperturbed HF response for the 
three considered parameter sets. The uniformity of this 
difference in the centroid energies suggests that there is 
only little change in the isoscalar residual interaction in 
the studied range of A's, and the positive shift due to the 
increase of the major shell gap is the dominant effect. 

The ISM response for Vsrg+DDI with A = 1.78 fm _1 
and Gogny D1S is rather similar. This can be seen as con- 
firmation of the effective low-momentum nature of the 
Gogny functionals. Furthermore, the similarity of the 
response also shows that there is a decoupling between 
the static and collective dynamical properties of the nu- 







mi /m [MeV] 


A[fm- X ] 


C 3 iv(A)[GeVfm 6 ] 


HF QRPA 


1.78 


4.95 


23.5 15.2 


2.02 


4.35 


22.5 14.2 


2.40 


3.42 


21.3 12.9 



TABLE III. Centroid energies of the isoscalar monopole 
strength distribution in 120 Sn for different Vsrg+DDI. The 
energy integration intervals for m /i were 5 — 50 MeV for the 
unperturbed HF response and 5 — 25 MeV for QRPA (see 
text). 



cleus, because the HFB ground-state energies obtained 
with Gogny D1S are very close to experimental ground- 
state energies in the tin isotopes, whereas Vsrg+DDI 
underestimates the ground-state energies by several MeV 
per nucleon, strongly depending on the SRG parameter 
A. The ground-state energy is an absolute quantity, while 
excitations are relative quantities that primarily depend 
on energy differences, and consequently, a simultaneous 
global shift of the levels involved in any single-particle 
transition will not alter the excitation energy. Such a 
shift seems to account for the bulk of the ground-state 
energy differences between Vsrg+DDI and Gogny D1S. 

The present findings are in concordance with previous 
RPA studies of closed-shell nuclei with UCOM interac- 
tions. In Rcf. [11], in particular, the RPA was formulated 
using the exact RPA ground state, and it was shown that 
the RPA ground-state correlations have little impact on 
the response. The RPA ground-state energy, however, 
contains many-body corrections [47], including second- 
order MBPT diagrams that recover much of the missing 
ground-state energy due to the perturbative nature of 
SRG and UCOM interactions [6, 32]. These observations 
are further indications that higher-order many-body cor- 
rections to the ground state are indeed small if SRG or 
UCOM interactions are used, but a more detailed and 
quantitative study of this subject would be of great in- 
terest. 

We conclude this section by investigating the effect 
of the SRG evolution on the energy-weighted sum rules 
(EWSRs). For the isoscalar monopole and quadrupole 
channels, there are the well-known classical expressions 
(sec, e.g., Rcf. [13]) 

2h 2 p 2 

Sls (E0) = ^—(N(rl)+Z(r 2 p )) , (27) 
Sls (E2)= 2 -^(N(rl) + Z(rl)) , (28) 

where (r^/ p ) are the intrinsic neutron and proton mean- 
square radii of the HFB solution, and for the dipolc 
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FIG. 8. (Color online) Isoscalar ( ) and isovector monopole response ( ) of selected calcium isotopes for Vsrg+DDI 

with (A[fm _1 ], C3jv[GeVfm 6 ])= (2.02, 3.87). The discrete strength distributions have been folded with a Lorentzian of width 
r = 1 MeV. 



strength, we have the Thomas-Rcichc-Kuhn sum rule [13] 



S(El) = 



2m 47r A 



(29) 



Equations (27) to (29) are derived by assuming that 
the interaction commutes with the transition operators, 
which is only the case for local interactions without 
isospin exchange contributions. Thus, the deviation of 
the calculated energy-weighted sum from the classical 
expressions is a measure for the size of momentum- 
dependent and isospin exchange contributions in the em- 
ployed Vsrg+DDI. 

In Fig. 7, we show the running energy- weighted sum 
in the isoscalar monopole and quadrupolc (ISQ), and the 
isovector dipole (IVD) channels. As discussed in Sec. 
II B, the latter also corresponds to the running energy- 
weighted El strength at low momentum transfer. The 
running sums exhibit general features that reflect the 
findings for the ISM response function shown in Fig. 
6, i.e., as A is lowered, the strength becomes less frag- 
mented, causing fewer jumps in Fig. 7, and the running 
starts at a higher excitation energy. 

At excitation energies of ~ 30 MeV, the isoscalar sum 
rules saturate, and Eqs. (27) and (28) are almost com- 
pletely exhausted. In general, the quadrupolc response 
exhibits a much weaker A-dependence than the monopole 
and dipole response, which is presumably due to the sur- 
face excitation character of the former, while the latter 
are volume modes. 

The total IVD sum exhibits the most pronounced A- 
dependence of the three investigated cases. For A = 
2.40 fm" 1 , the dipole EWSR is enhanced by 60% over the 
TRK value obtained from Eq. (29), and since strength is 
shifted to momentum-dependent terms as A is lowered, 
the enhancement increases to 70% for A = 1.78 fm -1 . 



While enhancement factors vary strongly among the 
available Skyrme and Gogny functionals [1, 48-50], our 
less phenomenological approach consistently favors a lim- 
ited range of values close to the enhancement factors ex- 
tracted from experimental photoabsorption cross-section 
data, which arc 70-75% for A > 100 [51]. 



B. 0+ Channel 

Having established the general effects of the SRG evo- 
lution on our QRPA results in the previous section, we 
now discuss the response of the calcium isotopic chain to 
the transition operators defined in Sec. II B. Figure 8 
shows the isoscalar and isovector monopole response for 
selected calcium isotopes. Since the results for the three 
SRG parameters are qualitatively similar, we display only 
the response for Vsrg+DDI with A = 2.02 fm -1 and 
C 3N = 3.87 GeVfm 6 (cf. Tab. I). 

The response exhibits the well-defined isoscalar giant 
resonance, which is found in the energy interval from 
10 to 20 MeV in the light calcium isotopes, and shifted 
to slightly lower energies as N increases. In the isovec- 
tor channel, we observe a distinct giant resonance peak 
around 32 MeV in 40 Ca, which fragments into a broad 
distribution with growing neutron excess. Beyond the 
major shell closure in 40 Ca, pronounced peaks start to 
emerge in a region around 10 MeV. For the lowest peak, 
in particular, the isoscalar and isovector response match 
closely in all isotopes up to 60 Ca, indicating that the re- 
sponse is strongly dominated by neutron transitions. As 
an illustrative example, we show the proton and neutron 
transition densities of the three major isoscalar peaks in 
50 Ca in Fig. 9. For the state at E = 9.352 MeV, the pro- 
ton transition density 8p p is negligibly small compared to 
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FIG. 9. (Color online) Proton ( ) and neutron ( ) 

transition densities for the major monopole peaks in 50 Ca. 
The light gray dashed lines indicate the calculated point- 
nucleon root-mean-square radius. 



the neutron transition density Sp n , and 8p n extends far 
into the nuclear exterior. The state at E = 12.204 MeV 
is of transitory character between the neutron-dominated 
low-lying state and the almost completely isoscalar Giant 
Monopole Resonance state at E = 17.303, where Sp p and 
5p n have similar extensions. 

While the low-lying strength emerges naturally and be- 
comes enhanced as we move toward the neutron drip line, 
the existence of low- lying strength in 48 Ca in our calcu- 
lations is a result of the insufficient spin-orbit splitting 
of the neutron 0/ levels, which creates a major shell for 
20 < N < 40 and turns 48 Ca into an open-shell nucleus 
(cf. Sec. IV). The otherwise qualitatively similar results 
of the recent Skyrmc QRPA survey by Terasaki and En- 
gel [52] do not exhibit a low- lying peak in 48 Ca since the 
EDFs used in their approach produce the correct major 
shell closures at N = 20 and 28. Furthermore, we ob- 
serve an overall enhancement of the low-lying excitations 
compared to Skyrme QRPA, likely due to differences in 
the single-particle spectra. In the Skyrme QRPA cal- 
culations, roughly 10-15% of the the isoscalar monopole 
EWSR are exhausted by states below 10 MeV for iso- 
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FIG. 10. (Color online) ISM centroids (top) and energies of 
the first excited + states (bottom) in the calcium isotopic 
chain, for Vsrg+DDI with (A[fin _1 ], C 3 jv[GeVfm 8 ])= (2.40, 
2.94) (•), (2.02, 3.87) (■), and (1.78, 4.41) (♦). Moments of 
the strength distribution were calculated in the energy inter- 
val 5 — 40 MeV. Experimental data are indicated by black 
symbols: Centroids for the 40 Ca giant monopole resonance 
were taken from the two analyses in Ref. [53], 0+ energies 
fromRef. [54]. 



topes with N > 30 (cf. Fig. 5 of Ref. [52]), while we 
obtain a growing exhaustion between 10% in 50 Ca and 
20% in 60 Ca in our calculations. Keep in mind, however, 
that the details of the running energy-weighted sum de- 
pend on the SRG parameter A, as discussed in Sec. V A. 

The upper panel of Fig. 10 shows the centroids of the 
ISM response for the three different sets of A's and C3at's. 
The behavior of the centroid energies as a function of 
the neutron number is identical for all three Vsrg+DDI, 
only the energies themselves increase as A is lowered, as 
discussed in Sec. VA. The centroids clearly reflect the 
shift of the ISM strength to lower energies as the neutron 
excess grows. The trend from 36 Ca to 40 Ca indicates a 
similar shift in proton-rich calcium isotopes. 

In the lower panel of Fig. 10 we display the energies 
of the first excited + states. We use the criterion intro- 
duced by Terasaki ct al. in Ref. [57] and discard states 
for which 



AA 



ij 



> 2. 



(30) 



because these states have pair-transfer character, and 
carry negligible strength from particle-hole excitations 
due to the consistency of our QRPA framework. In the 
closed-shell nuclei 40 Ca and 60 Ca (which has a closed- 
shell structure in our calculations), we find no low- lying 
+ state. In the case of 40 Ca, experimental data sug- 
gest that this state is intrinsically deformed [58] and can 
therefore not be described by the present spherical QRPA 
approach. 
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and isovector dipole response ( ) of selected calcium isotopes for Vsrg+DDI 

The discrete strength distributions have been folded with a Lorentzian of width 



FIG. 11. (Color online) Isoscalar ( — 
with (Alfm" 1 ], C 3 n{ GeVfm 6 ])= (2.02, 3.87). 
r = 1 MeV, and the isovector strength distribution has been scaled by a factor 250 for visibility. The arrows indicate the main 
giant dipole resonance peaks from photo-absorption cross-section data [55, 56]. 



The degree of A-dependence of the + states varies with 
the considered isotopes, and appears to be connected to 
the collectivity of each state. For the isotopes with 26 < 
N < 36, the collectivity of the + states is enhanced 
by pairing correlations, and their energies depend only 
weakly on A in the studied parameter range. The weak 
variation of the calculated + energies with the neutron 
number is qualitatively compatible with the experimental 
+ energies in the open-shell calcium isotopes between 
the physical major shell closures at N = 20 and N = 28. 
Thus, we expect to reproduce the experimental trends 
with a realistic interaction which yields the correct spin- 
orbit splittings. 



C. 1 Channel 

In Fig. 11, we show the isoscalar and isovector dipole 
response of selected calcium isotopes for Vsrg+DDI 
with A = 2.02 fur 1 and C 3N = 3.87 GeVfm 6 . The 
isoscalar response has two major components that cor- 
respond to lhid and 3hu excitations, respectively. For 
N > 20, these components are well separated energeti- 
cally by the isovector giant dipole resonance (IVGDR), 
i.e., states lying between 15 and 20 MeV carry signif- 
icant isovector dipole strength, but only little isoscalar 
dipole strength. In 36 Ca, however, the two states at 
15.397 and 20.414 MeV carry both. Their proton and 
neutron transition densities are displayed in Fig. 12. We 
immediately note that the roles of protons and neutrons 
are roughly reversed in the two states, while the densi- 
ties of the states are structurally similar: One density 
exhibits a pronounced peak near the surface, while the 
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FIG. 12. (Color online) Proton ( ) and neutron ( ) 

dipole transition densities for selected states in 36 Ca (see 
text). The light gray dashed lines indicate the calculated 
point-nucleon root-mean-square radius. 



other is wide and fiat. Due to the mismatched shapes of 
5p n and 8p p in the surface region, both states contribute 
significantly to the isoscalar and isovector response. 

Returning to Fig. 11, we note that the isoscalar re- 
sponse is further split into two components in the two 
previously identified regions. A low-lying peak around 
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FIG. 13. (Color online) Centroids of the IVD response in cal- 
cium isotopes for Vsrg+DDI with (A[ fin" 1 ], C 3N [ GeVfm 6 ])= 
(2.40, 2.94) (•), (2.02, 3.87) (■), and (1.78, 4.41) (♦). Mo- 
ments rrii were calculated in the energy intervals 10 — 40 MeV 

( ) and — 40 MeV( ), respectively. Experimental 

GDR peaks from photoabsorption data [55, 56] are indicated 

by (+)• 
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FIG. 14. (Color online) Proton ( ) and neutron ( ) 

dipole transition densities for the low-lying state at E = 
7.787 MeV in 44 Ca. The light gray dashed line indicates the 
calculated point-nucleon root-mean square radius. 



10 MeV is almost completely produced by proton tran- 
sitions in 36 Ca and 38 Ca (not shown). In Ca, neutron 
transitions start to contribute about 30% of the state's 
norm, making this mode largely isoscalar [59] , while be- 
yond the N = 20 shell closure, neutrons completely domi- 
nate this lower portion of the low-lying strength. The sec- 
ond peak at ~ 13 MeV that emerges beyond N = 20 also 
contains proton contributions, starting around 30 — 35% 
in 42 Ca and decreasing below 10% in 60 Ca. The peaks 
grow more pronounced and eventually overlap with grow- 
ing neutron excess, and the centroid of their strength is 
shifted to lower energies. 

In the high-lying region from 20 to 35 MeV, there is 
a very broad structure corresponding to the isoscalar gi- 
ant dipole resonance (ISGDR). While Fig. 11 suggests 
that there are two distinct groups of peaks, the states of 
both groups contain balanced contributions from proton 
and neutron transitions, as expected for the ISGDR, and 
there are no indications of major structural differences. 

As already mentioned above, the isovector dipole re- 
sponse of all calcium isotopes exhibits the pronounced 
giant dipole resonance (GDR) between 15 and 20 MeV, 
and the energies of the main resonance peaks agree rea- 
sonably well with experimental values extracted from 
photoabsorption data [56] . The corresponding GDR cen- 
troid energies, which are shown in Fig. 13, vary only 
weakly with N. The centroids exhibit the A-dependence 
we expect from our discussion in Sec. VA, i.e., the 
centroid energy increases as A is lowered from 2.40 to 
1.78 fm" 1 . 

With increasing neutron excess, the IVD response 
develops a pronounced low-lying peak structure below 
10 MeV. A similar but less pronounced peak is observed 
in the proton-rich nuclei 36 Ca (cf. Fig. 11) and 38 Ca 
(not shown). The emergence of low- lying El strength, 
or pygmy dipole strength as it is usually called, is high- 
lighted in Fig. 13 by including the centroid energy for the 
total IVD strength up to 40 MeV rather than just the re- 
gion of the GDR. It is a robust feature under variations of 



the SRG parameter A, but the details may differ because 
strength is shifted beyond 10 MeV as A decreases. 

For the low-lying El strength of 40 Ca, 44 Ca, and 48 Ca, 
experimental data are available from measurements con- 
ducted at the S-DALINAC [60]. In general our calcu- 
lations overestimate the energy of the PDRs and tend 
to underestimate the strength below 10 MeV. The over- 
estimation of the states' energy is ultimately in line 
with the conclusions from other studies in the literature, 
namely that in order to properly describe the low-lying 
El strength one needs to go beyond (Q)RPA and include 
particle-phonon coupling, as in, e.g., Rcf. [60]. 

From 42 Ca to 60 Ca, states below 10 MeV carry be- 
tween 2% and 4% of the El EWSR, corresponding to 
4%-7% of the TRK sum rule. The percentage of the 
sum rule increases smoothly with growing neutron excess, 
as expected at the QRPA level [60]. Almost all of the 
low-lying states carry both isoscalar and isovector dipole 
strength. As indicated in the discussion above, 70-80% 
of the norm of the low-lying states in 40 Ca is made up of 
proton transitions, while the low-lying states of the nu- 
clei with N > 20 are completely dominated by neutrons, 
which contribute 95-99% of the individual states' norms. 
Further analysis also shows that typically 10-20 different 
transitions contribute 1% or more each to the norm of 
each low-lying state, which confirms their collective na- 
ture, and supports their interpretation as pygmy dipole 
modes. The proton and neutron transition densities of 
the strongest low-lying state in 44 Ca that are shown in 
Fig. 14 exhibit the pygmy mode characteristics: 5p p and 
Sp n are in phase at short ranges, and out of phase just 
beyond the surface region of the nucleus, where 5p n also 
has a pronounced neutron tail, in agreement with other 
studies [61-63]. 
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FIG. 15. (Color online) Isoscalar ( ) and isovector quadrupole response ( ) of selected calcium isotopes for Vsrg+DDI 

with (A[fm _1 ], C3iv[GeVfm 6 ])= (2.02, 3.87). The discrete strength distributions have been folded with a Lorentzian of width 
r = 1 MeV. The arrow indicates the experimental centroid from Ref. [53]. 



D. 2+ Channel 

In Fig. 15, we show the isoscalar and isovec- 
tor quadrupole response of selected calcium isotopes 
for Vsrg+DDI with A = 2.02 fm" 1 and C 3N = 
3.87 GeVfrn 6 . We find the pronounced isoscalar giant 
quadrupole resonance (ISGQR) between 15 and 20 MeV, 
which exhibits moderate fragmentation. The isovector 
response is very broad in comparison, with the bulk of 
the strength residing in the IVGQR above 25 MeV. 

In the upper panel of Fig. 16, we show the centroid 
energies of the ISGQR for the three studied SRG param- 
eters. While the centroids exhibit the same trends under 
variation of A as the ISM and IVD centroids, we note 
that the resulting change in energies is smaller, covering 
an interval of roughly 1 MeV compared to the ~ 2 MeV 
range in the other cases. This observation matches the 
reduced sensitivity of the ISQ EWSR to variations of A 
discussed in Sec. V A. We also note that contrary to the 
ISGMR and IVGDR cases, where we typically approach 
the experimental centroid or peak energies from below as 
A is lowered, we overestimate the experimental centroid 
energy of the ISGQR in 40 Ca already for A = 2.40 fin -1 . 

Figure 15 exhibits prominent peaks at low energies, 
corresponding to the lowest 2 + states of the respective 
nuclei. In the lower panel of Fig. 16 we show the ener- 
gies of these states, which have been distinguished from 
states with pair-transfer character using the criterion (30) 
(cf. Sec. VB). A recent survey using a spherical Skyrme 
QRPA approach can be found in [57] . Similar to the first 
excited + states, the lowest 2 + states are insensitive to 
variations in the SRG parameter A, but for the latter, 
this also holds in the isotopes next to the major shell 
closures (cf. Fig. 10). In 40 Ca and 60 Ca, we do not 
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FIG. 16. (Color online) ISQ centroids (top), and energies 
of the first 2 + states (bottom) in the calcium isotopic chain, 
for Vsrg+DDI with (Affm" 1 ], C 3 iv[GeVfm 6 ])= (2.40, 2.94) 
(•), (2.02, 3.87) (■), and (1.78, 4.41) (♦). Moments of the 
strength distribution were calculated in the energy interval 
8 — 30 MeV. Experimental data are indicated by black sym- 
bols: the centroid for the 40 Ca giant quadrupole resonance 
was taken from Ref. [53] and 2+ energies from Ref. [54] . The 
2 + assignment is still tentative for the level in 52 Ca. 



find 2 + states below 10 MeV. The absence of the exper- 
imentally observed level in 40 Ca suggests that this state 
has a structure that cannot be described properly by our 
spherical QRPA. 

We find reasonable agreement between our results and 
the experimental 2 + energies in the open-shell calcium 



15 



0.02 



7 -o.i 



-0.3 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 
_ f\ \ E = 11.376 MeV _ 
f : \ : 

' rx \ 


j V 


1 " ^ ^ 

\ 
\ 

\ / 

\ / 54 Ca - 

i , , ,W, , i , , , i , , , : 


i 


| i 1 1 1 | i 1 1 1 J_J_L J-J-U'_l J_ 1 - 


X — x 
s \ 
s \ 

- \\ 
\\ 

; \\ 
\ 

■ T , , , 1 , , , , 


/ / 

I ' 
1 

/ / 
/ / 

/ ' 

V / ' 

\y 1 E= 18.554 MeV - 

\ / 54 Ca : 
1 >/ - 



r -i "i 

2 4 6 8 10 12 

r [ fm] 



FIG. 17. (Color online) Proton ( ) and neutron ( ) 

quadrupole transition densities for a Pygmy (top) and a GQR 
mode (bottom) in 54 Ca. The light gray dashed lines indicate 
the calculated point-nucleon root-mean square radius. 



isotopes. In semimagic spherical open-shell nuclei, the 2 + 
states lie at nearly constant energies along isotopic chains 
due to pairing correlations [64]. In our calculations, all 
isotopes from 42 Ca to 58 Ca exhibit this characteristic be- 
havior because the insufficient spin-orbit splitting of the 
0/ levels creates a major shell with 20 < N < 40. Conse- 
quently, we do not reproduce the increased experimental 
2 + energies in 40 Ca and 48 Ca, which signal the natu- 
ral major shell closures, nor the experimentally observed 
level at E = 2.563 MeV in 52 Ca, which has been tenta- 
tively identified as a 2 + state and whose increased energy 
is interpreted as a sign of an enhanced sub-shell clo- 
sure [65]. 

We conclude our discussion of the quadrupole response 
by focusing on the peak near 10 MeV, which emerges 
in the quadrupole strength distributions of neutron-rich 
calcium isotopes (cf. Fig. 15). All of the discrete 
states between 8 and 12 MeV contributing to this peak 
carry isoscalar and isovector strength in almost equal 
amounts because they are dominated by neutron transi- 
tions. They are also collective, receiving significant con- 
tributions from about 10 different transitions each, sug- 
gesting that they can be interpreted as pygmy quadrupole 
modes, whose existence was recently proposed in Ref. 
[66]. 

In Fig. 17, we compare the transition densities of the 
strongest 2 + state from the peak region to those of the 
main ISGQR state in 54 Ca. The transition densities of 
the ISGQR are in phase up to ~ 7.5 fm, and have a simi- 
lar magnitude. Their cxtrema are located within 1 fm of 
the calculated point-nucleon radius of 54 Ca. The proton 
transition density of the pygmy mode has a very small 



amplitude but similar extension to that of the ISGQR. 
The neutron transition density, on the other hand, ex- 
tends much further, and the bulk of its contribution to 
the transition operator is generated at much larger dis- 
tances than for the ISGQR state, which supports its in- 
terpretation as a neutron-skin excitation. 



VI. CONCLUSIONS 

We have developed a fully consistent QRPA frame- 
work for arbitrary AW interactions, represented by their 
two-body matrix elements, as well as 3A contact inter- 
actions via equivalent linearly density-dependent inter- 
actions. We use an intrinsic Hamiltonian and treat the 
Coulomb interaction exactly. The QRPA is built on the 
HFB ground states obtained with the code described in 
Ref. [27]. Since we do not truncate the QRPA configura- 
tion space and use the same Hamiltonian in the HFB and 
QRPA calculations, we achieve an excellent decoupling of 
the spurious strength. 

In the present work, we have employed our QRPA 
framework to study the nuclear response using NN in- 
teractions derived from Argonne V18 by means of the 
Similarity Rcnormalization Group, supplemented by a 
density-dependent interaction whose strength C$n is fit 
to the charge radii of closed-shell nuclei. By refitting C^n 
for each value of the SRG parameter A, we have been able 
to absorb the effect of the SRG evolution on the charge 
radii at the HFB level. 

Although the HFB ground-state energy strongly de- 
pends on A, and binding energy on the order of 4 — 
6 MeV per nucleon is missing, the theoretical single- 
particle spectra around the Fermi surface are reasonable, 
aside from the underestimation of the spin-orbit split- 
tings, which indicates the need for additional spin-orbit 
strength from the 3N interaction. Experimental odd- 
even binding energy differences, the analogs of pairing 
gaps in finite nuclei, are reasonably close to experimental 
values. Likewise, we obtain reasonable properties for the 
collective dynamics of the studied nuclei. 

In the nuclear response, we have identified two basic 
effects of the SRG evolution: The response is shifted to 
slightly higher energies as A is lowered, plausibly because 
interaction strength is shifted to momentum-dependent 
terms in the interaction, resulting in a reduction of the 
effective mass, and an increased spreading of the ma- 
jor shells in the single-particle spectrum. At the same 
time, the levels within a major shell are bunched closer 
together, resulting in less fragmentation at the (Q)RPA 
level. The important energy-weighted sum rules are sat- 
isfied for the studied SRG interactions, and exhibit only 
weak dependence on the SRG parameter A. Interest- 
ingly, we found that the enhancement of the isovector 
dipole EWSR over the classical Thomas-Reiche-Kulm 
value seems to be consistently close to experimental data, 
whereas the theoretical values strongly depend on the 
used EDF in Skyrmc of Gogny (Q)RPA studies. 
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Our results for the monopole, dipole, and quadrupolc 
response of the calcium isotopes are comparable to ex- 
isting studies based on phcnomenological EDFs. We 
have obtained reasonable centroid energies for the impor- 
tant giant resonances and identified multiple low-lying 
states in the dipole and quadrupole response that ex- 
hibit the characteristics of Pygmy modes. In accordance 
with other studies in the literature, the energies of these 
states are overestimated at the (Q)RPA level, illustrating 
the need to include many-body effects like quasiparticlc- 
phonon coupling to achieve a proper description of low- 
lying strength. 

The next major development stage of our QRPA 
framework will be the switch to SRG-evolved chiral NN 
and 3iV interactions. Since the latter contain spin-orbit 
and tensor structures, we expect an impact on the issue 
of the underestimated spin-orbit splittings, although it 
remains to be seen whether the problem will be fixed. 
Realistic spin-orbit splittings are essential to allow stud- 
ies of heavier isotopic chains, where the ground states de- 
velop a deformation for the currently used Vsrg+DDI, 
causing the breakdown of the spherical QRPA. 

The extension of our HFB+QRPA framework to de- 
formed nuclei is an obvious direction for further research, 
but the treatment of realistic NN and 37V interactions 
in deformed bases will require significant work, particu- 
larly on the transformation between relative partial waves 
and uncoupled deformed single-particle states. Besides 
the treatment of deformation, we are looking into exten- 
sions toward more complicated configurations like 4qp ex- 
citations, i.e., a Second QRPA analogous to the Second 
RPA described in Ref. [12], and quasiparticle-phonon 
coupling. 

Our initial QRPA survey of the calcium isotopes has 
indicated interesting structural features like low-lying 
modes, including the prominent IHlj peak in the isoscalar 
dipole response, which was studied for closed-shell nuclei 
in Ref. [59] . More detailed studies of the physics of these 
modes in open-shell nuclei will be the subject of research 
in the near future. 
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Appendix A: Explicit Form of the QRPA Matrices 



The matrices A and B are given by [fj, = (n^l^j^r^)] 



1 



f- 0„„i l. 



+ F(^fi'vv'; J) (u^VpiUvVv' + (u «-» v)) 

- (-iy" +j "'- J F(nn'u'v; J) (u^v^V + ( u v )) 

+ G(fj,fi'vi/'; J) (u^u^UyUy, + (u «-» v)) j (Al) 



and 



x ji^/i/iW; J) [y^u^UyV^ + (u o v)) 

- (-lji'+W-^^Vi/; J) (UftVft'UvVv, + (ll <-» v)) 

— G{[ip!vv '; J) {u^u^VvVv' +{u^v))X. (A2) 
The single-quasiparticle term Hj^, reads 

(A3) 

where A is the chemical potential (separate for neutrons 
and protons), 

K»> = fi - jj W + + !)w^ ( A4 ) 



is the particle-hole field, and 



(JbfJb' VV U'V U V 



(A5) 



is the pairing field in the canonical basis representation. 
Note that both h and A arc diagonal in I and j and inde- 
pendent of magnetic quantum numbers due to parity and 
spherical symmetry, and the two-body matrix elements 
V/ifi'vu' contain all two-body terms of the Hamiltonian, 
i.e., nuclear and electromagnetic interactions, as well as 
the two-body center-of-mass correction [16, 27]. 

We adopt the notation of Ref. [18] for the particle- hole 
and particle-particle channel matrix elements, which are 
given by 

F(wiW; J) = £(-l)v+^'( 2 j' + i) ( -V J \ 

(A6) 



and 



x (fjLj/;J'\v\n'v;J') 



G(/i/iVi/; J) = (fJ,fjf] J\ v \vv'\ J) 



(A7) 



As stated in Sec. II A, the two-body matrix elements are 
antisymmetrized but not normalized. 
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Appendix B: Density-Dependent Interaction & 
Rearrangement Terms 

The contribution of the density-dependent interaction 
(23) to the energy is given by 

_ _ (*\v{ P }\y) _ (f|Ei<,-%-[p]l*) 



C 3N n / drr 2 p(r)p n (r)p p (r) 



(Bl) 



where 



fh(r) = 



2j ^±vl\R»{r)\\ (B2) 



Rfj,{r) are the canonical wavefunctions, and the total den- 
sity is 



P(r) = Pn(r) +p P (r). 



(B3) 



Note that there is no contribution to the pairing energy, 
since the interaction acts only in the (S, T) = (1, 0) chan- 
nel, i.e., like-particle matrix elements vanish. 

Variation of Eq. (Bl) w.r.t. to the density matrix 
yields the particle-hole field that has to be added to h^i 
[Eq. (A4)]: 



r 



4 



drr 2 R;(r)p(r) [p(r) - p T »] R^(r) 



' drr 2 R*(r)p p (r)p n (r)R^(r), (B4) 



where the second line is the rearrangement term. 

In the QRPA matrices, the DDI contributes the matrix 
elements 



{pm^vm v \ v ph [p] Ip'm^v'rriui 



2 E o 



(B5) 

to the particle-hole interaction [33] , which can be plugged 
into Eq. (A6). Denoting the normal and rearrangement 
matrix elements Fq and Fi , respectively, one obtains for 
identical isospins (r = t m = 7y = t„ = r„< ) 



F o (pp J 'vv';J) = 0, 



(B6) 



and (j = V27TT) 

X (-l)^-J, (1 + (1 + 

x J dr r 2 ^ (r)R„> (r)R v (r)R v > (r) |p(r) - p T (r) . 

(B7) 

If the isospins are only pairwise identical (r^ = Ty ^ 
t v = t v >), 

F (pp'i/i/'; J) 



12 



K(pp!vv' ; J) (-1)^ +J« 



{(i4v-IIJ°>0'4>'-il J o) 

x (_!)v-J> ((-1) J+WI '' +2) 
+ (i4v|l J1 )0'4>'llJ 1 >(- 1 ) V+i '"} . (B8) 



and 



= % L oviv-iiw4v-||Jo) 

x (_i)JW» (l + (_i)'m+V+ j ) (l + (_i)'-+'^+J) 
x K(pp!w'\ J) , (B9) 

where the radial integral is given by 

K{ppvv\ J) 



AttJ 2 



JoJW I drr 2 R ^ r)R ^ {r ) Ry (r)R ul ( r )p(r) 



(BIO) 



We have confirmed that these matrix elements match the 
more general versions derived in [33] analytically and nu- 
merically. 
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